Discovery of VHE gamma-ray emission from 1ES 121 8+30.4 
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ABSTRACT 

£3 ■ 

, The MAGIC collaboration has studied the high peaked BL-Lac object 1ES 1218+30.4 at a redshift z = 0.182, using the MAGIC imag- 
• • o ing air Cherenkov telescope located on the Canary island of La Palma. A gamma-ray signal was observed with 6.4 cr significance. The 
. m differential energy spectrum for an energy threshold of 120 GeV can be fitted by a simple power law yielding F E (E) = (8.1 ± 2.1) • 
, 10~ 7 (is/250 GeV)~ 3 - 0±0 - 4 TeV~' irr 2 s -1 . During the six days of observation in January 2005 no time variability on time scales of days was 
■ found within the statistical errors. The observed integral flux above 350 GeV is nearly a factor two below the the upper limit reported by the 
d ; Whipple Collaboration in 2003. 
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1. Introduction 

Since the detection of blazars as sources of high energy 
gamma-rays by EGRET on board the CGRO, the search for 
the very high energy (VHE) counterparts of these sources 
has been a major goal for ground-based gamma-ray astron- 
omy. The first blazar detected at VHE energies was Mrk421 
JPunch et al. 19 92 1. Subsequent to the discovery of Mrk421, 
only a small number of extragalactic VHE sources was found, 
although the sensitivity of ground-based telescopes is superior 
in their energy range to that of EGRET. The total number of 
VHE blazars reported in the literature as of February 2006 
amounts to 11, while the number of high-confidence blazars 
in the 3rd EGRET catalogue is 66 ( IHartman et al. 19991 1. 

Stecker et al. (1992) pointed out that the attenuation of 
gamma-rays due to photon-photon interactions with low en- 
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ergy photons from the extragalactic background radiation was a 
likely explanation of this deficit. In fact, the redshifts of blazars 
detected so far above 100 GeV are at rather low values as ex- 
pected from predictions of the correlation between the gamma- 
ray attenuation and the redshift of the source, known as Fazio- 
Stecker relation (IFazio & Stecker 19791 IKneiske et al. 2004> . 
These primary gamma-rays, absorbed by pair production, turn 
into secondary gamma-rays conttibuting to the extragalactic 
background light ( Aharo nian et al. 1994 1. 

The VHE gamma-ray emitting blazars known so far all be- 
long to the class of high-frequency -peaked BL Lacertae objects 
(HBLs, Padovani & Giommi 'l995> . a subclass of blazars char- 
acterised by a low luminosity when compared with quasars and 
a synchrotron peak in the X-ray band. Their Spectral Energy 
Distribution (SED) is characterized by a second peak at very 
high gamma-ray energies. In synchrotron-self-Compton (SSC) 
models, it is assumed that the observed gamma-ray peak is due 
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to the inverse-Compton emission from the accelerated elec- 
trons up-scattering previously produced synchrotron photons to 
high energies ( Maraschi et al. 199 2 1. In hadronic models, sec- 
ondary electrons and positrons, arising from photo-meson pro- 
duction, initiate electromagnetic cascades (Mannheim 1993 
IMuecke & Protheroe 200 1> . A compilation of blazars with 
known X-ray spectra allowing their classification as HBLs is 
given by |Donatoet al. (200 1)| 

During the first year of operation of the Major Atmospheric 
Gamma-ray Imaging Cherenkov (MAGIC) telescope, a sam- 
ple of X-ray bright (F\ k e v > 2 /Jy), northern HBLs at mod- 
erate redshifts (z<0.3), including 1ES 1218+30.4, which is 
located at a redshift z = 0.182 fVeron-Cetty & Ver on 2003 1 
have been observed. Predictions based on models invoking a 
synchrotron self-Compton ( Costamante & Ghisellini 2002 1 or 
hadronic ( Mannhei m et al. 19 96 1 origin of the gamma-rays and 
the shape of the time-averaged SED show that these HBLs 
should be detectable by imaging air Cherenkov telescopes 
(IACTs), although most of them are too faint in the EGRET 
energy range to were detected during its all-sky survey. 

For 1ES 1218+30.4, the EGRET upper limit is 
F( > 100 MeV) « 10 _1 1 erg cm~ 2 s . Several observations with 
the Whipple telescope between 1995 and 2000 resulted in an 
upper flux limit of F(> 350 GeV) = 8.3 • 10~ 8 photons nr 2 s" 1 
(corresponding to ~8% of the Crab Nebula flux) 
JHoran et al. 2 004 1. The source was also observed by 
HEGRA between 1996 and 2002. An upper flux limit above 
840 GeV of 2.67 • 10' 8 ph otons m" 2 s" 1 (or 12% o f the Crab 
Nebula flux) is reported by Aharonian et al. (2004) 

In this letter, we report the first detection of VHE gamma- 
rays (> 100 GeV) from the direction of 1ES 1218+30.4. In 
Sect. 2 we briefly discuss the observations and the data set 
while we describe in Sect. 3 the data analysis method and the 
results. In Sect. 4 we discuss the results in the context of the 
spectral energy distribution of 1ES 1218+30.4 across a broad 
energy range. 

2. Observations 

The MAGIC telescope is a single dish IACT, located on the 
Canary island of La Palma (28.8° N, 17.8° W, 2200 m a.s.l.). A 
17 m diameter tessellated parabolic mirror with a total surface 
of 234 m 2 , mounted on a light-weight space frame made from 
carbon fiber reinforced plastic tubes, focuses Cherenkov light 
from air showers, initiated by gamma-rays or charged cosmic 
rays, onto a 576-pixel photomultiplier (PMT) camera with a 
field-of-view of 3.5°. The analogue signals are transported via 
optical fibers to the trigger electronics and each channel is read 
out by a 300 MHz FADC. Further details on the telescope can 
be found in |Baixeras et al. (2004)| and |Cortina et al. (2005)| 

MAGIC observed 1ES 1218+30.4 from January 9th to 15th 
2005 at zenith angles (ZA) between 1.5° and 13° during six 
moonless nights, for a total observation time of 8.2 h. To de- 
termine the background, so-called off-data were taken in addi- 
tion, by pointing the telescope to a nearby sky region where no 
gamma-ray source was expected. The off-data cover the same 
ZA range with a similar night-sky background light intensity. 
For the analysis, 6.5 h off-data, taken between January 9th and 



January 11th, were used, which match the observation con- 
ditions and detector performance of the on-data. At the same 
time the KVA telescope (http://tur3.tur.iac.es/) on La Palma ob- 
served 1ES 1218+30.4 in the optical range. 

The signal was discovered with our automatic analysis 
JDorner et al. 20 05 ). The final analysis is decribed hereafter. 

3. Data analysis and results 

The data were processed using the MAGIC Analysis 
and Reconstruction Software (MARS) dBretz 2005at . A 
description of the different analysis steps can be found 
in |Gaug et al. (2005)| (including the calibration) and 
Bretz (2005b) The moments up to third order of the light 
distribution are used to characterise each event by a set of 
image parameters ( Hillas 1985). 

For background suppression, a SIZE-dependent parabolic 
cut in WIDTH x LENGTH is applied ( Riegel 2005). 
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Fig. 1. Sky map of the region around the position of 
1ES 1218+30.4 (cross) for gamma events with an en- 
ergy threshold of ~ 140 GeV. The scale is in units of 
10 events/arcmin 2 . 

To reconstruct the origin of the shower in the camera 
plane, the DISP method is employed (Lessard et al. 200T1 to 
estimate the distance between the centre of gravity of the 
shower and its origin. The third moment determines the di- 
rection of the shower development. The constant coefficient 
^ from the parameterisation of DISP in the original approach 
is replaced by £ + ft • (LEAKAGE) 6 , LEAKAGE being 
the fraction of light contained in the outermost camera pix- 
els. Thereby the truncation of shower images at the camera 
border is taken into account. These coefficients were deter- 
mined using simulated data. Simulated gamma-showers (MC) 
were produced by COR SIKA, version 6.023 JHeck et al. 19981 
Majumd ar et al. 200 5 1 for ZA below 20° and energies between 
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lOGeV and 30TeV, following a power law with a spectral in- 
dex -2.6. The cut coefficients for the background suppression 
were optimised using data of the Crab Nebula, taken at similar 
ZA in a time slot of several weeks around the observation. 

Figure^ shows the background-subtracted distribution 
of reconstructed shower origins centred at the position of 
1ES 1218+30.4. Each shower origin is folded with a two- 
dimensional Gaussian with a standard deviation of 0.06°, cor- 
responding to half of the sigma of the point-spread function ex- 
pected for a point-like gamma-ray source as seen by MAGIC. 
The centre of gravity of the observed excess coincides within 
the systematic pointing uncertainty of 0.04° with the nominal 
source position. 




Ra: 12h21m22 
Dec: 30°10'38 

Observation: 8.2h 
Significance: 6.4a 
Events: 

Excess: 560 

Background: 3308 

Off-Scale: 1.2 
Zenith Angle: 1°-13° 



Fig. 2. Distribution of 6 2 for on- and off-data (grey shaded). The 
signal region is marked by the dashed line. 
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Fig. 3. Differential energy spectrum of 1ES 121 8+30.4. The up- 
per limits correspond to a 90% confidence level. 



(grey-shaded region). The total systematic error for the spec- 
tral slope is estimated to be + 0.7 - 0.4. Additional systematics 
for the flux level are estimated to be in the order of ~40%. The 
main contributions are uncertainties of the atmospheric con- 
ditions, the mirror reflectivity and the effective quantum effi- 
ciency of the PMTs. 



Figure|2 shows the distribution of the squared angular dis- 
tance, 9 2 , between the reconstructed and the nominal source 
position, for both the on- and off-data. The background distri- 
bution is not flat due to the limited trigger area, which means 
that the highest acceptance of the detector is in the center and 
decreases outwards. Due to different observation times, the off- 
data were scaled to match the on-data in the region 0.14 < 
(6/°) 2 < 0.64, where no bias from the source is expected. The 
observed excess of 560 events has a statistical significance of 
6.4 standard deviations (according to lLi & Ma 1983 1. The en- 
ergy threshold, defined as the peak energy of simulated gamma- 
rays with a differential spectrum oc E~ 3 ' surviving all cuts, is 
-140 GeV. 

For the energy estimation, the Random Forest regression 
method was applied ( Breiman 200 1 1, trained by the above men- 
tioned MC sample, yielding an average energy resolution of 
24%. Figure|3]and table^show the reconstructed energy spec- 
trum. To avoid systematic effects of cut-efficiencies to the MC- 
sample, the acceptance window for gamma-rays is enlarged, 
yielding an energy threshold of ~ 120 GeV. Spill-over effects 
to neighbouring bins were corrected by the number ratio be- 
tween the energy distribution of a MC sample compared to the 
distribution of estimated energy. A power law is fitted to the 
measured spectral points yielding (x 2 INDF = 1.1/4): 

, E ,-3.0±0.4 

F E (E) = (8.1 + 2.1) • 10~ 7 rrT 2 s _I TeV^(l) 

V ; V ; \250GeV/ V ; 

The error bars (lcr) show the statistical uncertainty only. 
The systematic error resulting from the analysis was investi- 
gated by changing cut-coefficients and the initial MC spectrum 



Table 1. Differential flux of 1ES 1218+30.4 and the number of 
measured excess events together with the statistical error. For 
the flux also the systematic error is given (see text for details). 



E [GeV] 


F E (E) ± stat +sys -sys [m 1 


s-'TeV-'] 


Nexc ± stat 


87 


(1.5 + 1.2 + 8.7-1.2) 


io- 5 


27 + 21 


130 


(9.2 + 4.2 + 3.1-7.5) 


IO" 6 


163 ± 73 


190 


(1.8 + 0.9 + 0.8-0.4) 


10~ 6 


146 + 71 


280 


(4.3 + 3.1+0.9-1.4) 


io- 7 


64 + 46 


420 


(2.1 ± 1.1+0.4-1.1) 


io- 7 


61+30 


630 


(4.5 + 3.5 + 1.9-2.7) 


io- 8 


22 + 17 



The data were also analysed with two other, inde- 
pendent analysis techniques, using dynamical super-cuts 
(Kranich 1997) or the Random Forest method for background 
suppression, as used in previous MAGIC observations (e.g. 
Alb ert etal. 20 06 1. Within the statistical errors, the same sig- 
nificance, flux and differential spectrum were obtained. 

For the analysis of time variability, the sample was di- 
vided into six sub-samples, one for every night of observa- 
tion, with ~85 min observation time each. Figure^ shows the 
integral flux above 100 GeV for each night. Within the sta- 
tistical errors, we get a constant flux (x 2 /NDF - 2.4/5) of 
F(> 100 GeV) = (8.7 + 1.4) • 10~ 7 photons irT 2 s" 1 . The sys- 
tematic error is in the order of ~40% as discussed above. The 
integral flux is ~30% higher than one would expect from inte- 
gration of Eq. Q. This can be understood, taking into account 
that the flux at 130 GeV (peak of the event distribution) lies 
~35% above the fit value. 
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Fig. 4. Integral flux ( > 100 GeV) of 1ES 1218+30.4 for every 
sub-sample. The error bars (lcr) show the statistical uncertainty 
only. 



4. Discussion and conclusions 

The discovery of VHE gamma-ray emission from 
1ES 1218+30.4 provides further evidence that HBLs gen- 
erally exhibit gamma-ray emission at a luminosity comparable 
to the X-ray luminosity. Archival X-ray data show that the 
X-ray flux at 1 keV varies in the range 0.8 to 4.5 in units of 
10~ n ergcirr 2 s _1 (ROSATi, which is below the sensitivity 
of the All Sky Monitor (ASM) on board the RXTE X-ray 
satellite. The observed gamma-ray flux at 100 GeV is with 
(1.9+ 1.1) -10- 11 ergcirr 2 s" 1 at the same level as published 
X-ray fluxes. Contemporaneous observations from Swift in 
October 2005 showed a flux comparable to that measured by 
BeppoSAX in 1999 JGiommi 2006> . which is also within the 
range of the one measured by ROSAT. 

There are no indications from simultaneous optical obser- 
vations (KVA) of a flare to have occurred during the time of 
the gamma-ray observations, neither does the gamma-ray light- 
curve show signs of significant variability. The spectrum is con- 
sistent with upper limits at higher energies determined in the 
past. 

Taken the upper limit of EGRET into account, the gamma- 
ray peak of the SED is constrained to lie in the 1 GeV to 
100 GeV regime, unless attenuation of the gamma-rays by pair 
production in the metagalactic radiation field masks an emit- 
ted bump in the spectrum at higher energies. Detailed mod- 
elling of the attenuation thus becomes the task to unfold the 
emitted spectrum from the observed one. Sizeable attenuation 
is expected from current models of the evolving extragalac- 
tic (=metagalactic) background light dHauser & Dwek 20011 
IKneiske et al. 20041 . 

At hard X-ray energies, SSC models predict a sharp roll- 
over in the SED, in contrast to the predictions of hadronic 
models involving electromagnetic cascades. Observations in 
the hard X-ray domain, e.g. using Swift or Suzaku (Astro-E2), 
and in the intermediate gamma-ray regime with GLAST, will 
be crucial to provide information about the missing link be- 
tween the synchrotron peak and the gamma-ray peak. 
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